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, Abstract. Lichnerowicz's algebra of differential geometric operators acting on symmetric 

^ ■ tensors can be obtained from generalized geodesic motion of an observer carrying a complex 

[ tangent vector. This relation is based upon quantizing the classical evolution equations, and 

■ identifying wavefunctions with sections of the symmetric tensor bundle and Nocthcr charges 
with geometric operators. In general curved spaces these operators obey a deformation of the 
Fourier Jacobi Lie algebra of sp(2,M.). These results have already been generalized by the 
authors to arbitrary tensor and spinor bundles using supcrsymmetric quantum mechanical 
models and have also been applied to the theory of higher spin particles. These Proceedings 
review these results in their simplest, symmetric tensor setting. New results on a novel and 

■ extremely useful reformulation of the rank 2 deformation of the Fourier-Jacobi Lie algebra in 
fH | terms of an associative algebra are also presented. This new algebra was originally motivated 

by studies of operator orderings in enveloping algebras. It provides a new method that is 
superior in many respects to common techniques such as Weyl or normal ordering. 
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■ 1 Introduction 

■ The study of geometry using first quantized particle models has a long history. Notable examples 
C^- , are the study of Pontryagin classes and Morse theory in terms of jV = 1 and M = 2 supersym- 

metric quantum mechanical models [TJ [2]. The supercharges of those models correspond to 
Dirac, and exterior derivative and codifferential operators acting on spinors and forms, respec- 
tively. The model we concentrate on here describes gradient and divergence operators acting 
on symmetric tensors and therefore involves no supersymmetries at all. Hence, even though 
the symmetries of this model are analogous to supersymmetries, no knowledge of superalge- 
bras is required to read these Proceedings. All the above models fit into a very general class 
of orthosymplectic spinning particle theories studied in detail by the authors in [3]. Spinors, 
differential forms, multiforms 01 E El El [8] and symmetric tensor^ HE] are all fitted into 
a single framework in that work. Here we focus on the symmetric tensor case, both for its 
simplicity, and because we want to present new results on the symmetric tensor Lichnerowicz 
algebra developed in [TO] . 

The underlying classical system is geodesic motion on a Riemannian manifold along with 
parallel transport of a complex tangent vector. This is described by a pair of ordinary differen- 
tial equations to which we add further curvature couplings designed to maximize the set of 
constants of the motion. Of particular interest are symmetries interchanging the vector tangent 



*This paper is a contribution to the Proceedings of the 2007 Midwest Geometry Conference in honor of 
Thomas P. Branson. The full collection is available at |http://www.emis.de/journals/SIGMA/MGC2007.html| 

x See [TT] for a flat space discussion of the symmetric tensor theory and |12l I13| for its origins in higher spin 
theories. 
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to the manifold with the tangent vector to the geodesic. These are analogous to supersymmet- 
ries and correspond to gradient and divergence operators. This correspondence is achieved by 
quantizing the model. The complex tangent vector describes spinning degrees of freedom so that 
wavefunctions are sections of the symmetric tensor bundle. The Noether charges of the theory 
become operators on these sections. In particular, the Hamiltonian is a curvature modified 
Laplace operator. In fact, it is precisely the wave operator acting on symmetric tensors intro- 
duced some time ago by Lichnerowicz [9] on the basis of its algebraic properties on symmetric 
spaces |10j . Moreover, the set of all Noether charges obey a deformation of the Fourier-Jacobi 
Lie algebra sp(2,W) J . The classical model is described in Section [2j while its quantization and 
relation to geometry are given in Section [3l 

Applications, such as higher spin theories \10\ I14j. call for expressions in the universal en- 
veloping algebra U{sp{2, R) J ) involving arbitrarily high powers of the generators. Manipulating 
these expressions requires a standard ordering, oft used examples being Weyl ordering (aver- 
aging over operator orderings) or normal ordering (based on a choice of polarization such that 
certain operators are moved preferentially to the right, say). In a study of partially massless 
higher spins [15], we found a new operator ordering scheme to be particularly advantageous [10] . 
The key idea is to rewrite generators of the sp(2, R) subalgebra, wherever possible, as powers 
of Cartan elements or the quadratic Casimir operator. Immediately, this scheme runs into 
a difficulty, namely that the remaining sp(2, R) generators do not have a simple commutation 
relation with the quadratic sp(2, R) Casimir. This problem is solved by a trick: we introduce 
a certain square root of the quadratic Casimir whose role is to measure how far states are from 
being highest weight. Then we use this square root operator to construct modified versions of 
the remaining sp(2,R) J generators. Instead of a simple Lie algebra, we then obtain an elegant 
associative algebra, which we denote U(sp(2, R)), with relations allowing elements to be easi- 
ly reordered. This algebra is described and derived in detail in Section [H The final Section 
discusses applications and our conclusions. 



2 The classical model 

Let {M^q^jj) be an n-dimensional (pseudo-)Riemannian manifold with an orthonormal frame e m 
so thatj ds 2 = dx^g^ v dx u = e m rj mn e n . We consider the motion of an ant x^(t) - as depicted 
in Fig. Q] - who carries a complex vector z m (t) (expressed relative to the orthonormal frame) 
tangent to M. (In physics nomenclature, z m is referred to as commuting spinning degrees of 
freedom.) The Levi-Civita connection will be denoted by V. The ant determines its path and 
in which direction to hold the complex vector by the system of generalized geodesic ODEs 

J__ - m * n , \jHr>m r * n „* s 

Vz m 

= iR™ n r sZ n Z* r Z S . (1) 

dt 

The non-linear couplings to the curvature tensor on the right hand side of these equations have 
been carefully chosen to maximize the set of constants of the motion. They may obtained by 
extremizing a generalized energy integral 

f fl Vz m 1 1 

O — / (IIS —X g^i/X -+- lZ m — h — It n s Z m Z Z r Z > . yZ) 



2 Although the metric signature impacts the unitarity of the quantum Hilbert space of our model, all the results 
presented here hold for arbitrary signature. Similarly, none of our results depend on the existence of a global 
orthonormal frame. 
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To study constants of the motion, we look for symmetries of this action principle. The most 
obvious of these are translations of t — > t + £ along the parameterized path traversed by our ant. 
Infinitesimally this yields the invariance 

5x» = £fc", 5z m = iz m . 

Less trivial, are symplectic transformations of (z m ,zj^), 

5z m = -az m + (3z* m , 6z* m = 1Z m + az* m . 

The parameters (a,/3, 7) are real and correspond to the Lie algebra sp(2, R). The astute reader 
will observe that the symmetry transformation of z* m is not the complex conjugate of z m . 
Nonetheless, treating z rn and z* rn as independent variables, the above sp(2,M.) transformations 
do leave the action invariant. This is in fact sufficient to ensure existence of corresponding 
constants of the motion and Noether charges. In the quantum theory, these charges will play 
an important geometric role. 

The most interesting symmetries of the model interchange the complex tangent vector z rn 
with the tangent vector to the ant's path 

5x^ = i{z*^e - z^e*), Vz p = x»e. (3) 

Here, T> is the covariant variation and is defined by Dv 11 = 5v^ + Tp cr 5x p v' 7 where T denotes the 
Christoffel symbols. It saves one from having to vary covariantly constant quantities. 

The transformations ([3|) are not an exact symmetry for an arbitrary Riemannian manifold. 
In fact, the action ((2j) is invariant only when the locally symmetric space condition 

V«JVv* = °> ( 4 ) 

holds. Or in other words, the Riemann tensor is covariantly constant. Constant curvature spaces 
provide an, but by no means the only, example of such a manifold. 

To compute constants of the motion we work in a first order formulation 



This and the evolution equations ([TJ also follow from an action principle 

= Jdt L^x" + iz* m z m - \it^k u + \R m n \ z* m z n z;.z s } , (5) 
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where the covariant and canonical momenta ir^ and p^ are related by 

Here the spin connection is determined by requiring covariant constancy of the orthonormal 
frame V p e M m = d p e p m - T^e™ + Lo p m n e p n = 0. 

From the first order action ([5]) we immediately read off the contact one-form p p dx^ + iz^ n dz 
which is already in Darboux coordinates, so Poisson brackets follow immediately 

The Noether charges for the symmetries of the model can now be computed 

H — — TT p g^ 7T U 2^ n S ^ rn ^ i 



IZ^TT 



(6) 



The first of these is the Hamiltonian. We have arranged the symplectic symmetry charges in 
a symmetric matrix / using the isomorphism between the symplectic Lie algebra and symmetric 
matrices. The remaining charges appear as a column vector v since they in fact form a doublet 
representation of sp(2,K). It is important to remember that this latter pair of charges are 
constants of the motion for locally symmetric spaces only. 



3 Quantization and geometry 

Quantization proceeds along usual lines replacing the Poisson brackets by quantum commutators 
[p, x] = —ih and [z, z*\ = K. We set H = 1 in what follows and represent the canonical momentum 
as a derivative acting on wavefunctions ^(x^) 

1 d 



i dx^ 

The spinning degrees of freedom become oscillators acting on a Fock space. Rather than using 
the standard notation z m = a m annihilating a Fock vacuum a m \0), we represent |0) = 1 and to 
preempt their geometric interpretation, set 

d 

z*^ = dx 1 *, z 



" d{dx») 
Therefore, wavefunctions become 



s=0 



or in words - sections of the symmetric tensor bundle SM over M. Therefore, we can now start 
relating quantum mechanical operations to differential geometry ones on symmetric tensors. 
Firstly, the quantum mechanical inner product yields the natural inner product for symmetric 
tensors 



oo 



The Symmetric Tensor Lichnerowicz Algebra 



5 



Furthermore, the covariant momentum corresponds to the covariant derivative iir^ = V^VP 
(it is necessary to contract the open index fx with dx^ for this to hold true for subsequent 
applications of 7r M ). 

Next we turn to the symplectic symmetries / in ©. We call the off-diagonal charge 
d 

N = dx^- 



d{dx») ' 

which simply counts the number of indices of a symmetric tensor 

N Wi-ms^ 1 ■■■dx fl3 = sip lil ... fls dx^ ■■■dx^. 
We call the diagonal charges 

E = dx^dx\ ^=Q^/ W Q^ry 

as they produce new symmetric tensors by either multiplying by the metric tensor and sym- 
metrizing, or tracing a pair of indices 

g ^m.-.^dx^ ■ --dx^ = g illl ^ lls ... lls+2) dx f11 ■ --dx^ 2 , 

tr • • • dx>ls = s ( s ~ l^Vi-Ms-a^ 1 • • • dx^-\ 

These three operators obey the sp(2,R) Lie algebra 

[N,tr] = -2tr, [N,g]=2g, [tr, g] = 4 N + 2n. 
We call its quadratic Casimir 

c = gtr-N(N + ?i-2). 

Trace-free symmetric tensors with a definite number of indices, (N — s)^f = = tr^, are the 
highest weight vectors for unitary discrete series representations of this sp{2, K) algebra. 

The Noether charges v in ([6]) are linear in momenta and therefore covariant derivatives, when 
acting on wavefunctions. We call them the gradient and divergence, 

grad = div = ^ V, 

because they are natural generalizations to symmetric tensors of the exterior derivative and 
codifferential for differential forms. To be sure 

grad Wi...^^ 1 • • • dx^ = V (w Vv2...^ + i)^ w ■ ■ ■ dx^+\ 
div^...^^ 1 • • • dx"' = sV^ 1 ... (1 ^ 1 d/ 1 • • • da?- 1 . 

This pair of operators forms the defining representation of sp(2, R) 

[N,grad] = grad, [N, div] = —div, [tr, grad] = 2 div, [div,g] = 2 grad. 
It remains to commute the operators div and grad. The result is 

[div, grad] = A - R** , (7) 
where A = V^V^ is the Bochner Laplacian and 

R** = ^p^^Q^^d^y 
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This relation is closely analogous to that for the exterior derivative and codifferential {d, 5} = 
Ap where Ap is the form-Laplacian. Here, since we are dealing with symmetric tensors, the 
anticommutator is replaced by a commutator. A shrewd reader might sense that the second order 
operator on the right hand side of ([7]) should be related to the quantum mechanical Hamiltonian 
operator. This is indeed the case; calling □ = A + R^^, we have 

[div,grad] = □ - 2R** , 

where □ = — 2H so long as an appropriate operator ordering is chosen for the Hamiltonian 
(a full account is given in [3]). On any manifold 

[□,g] = [D,N] = [D,tr] =0. 

Moreover, whenever the symmetric space condition (J4|) holds, the operator □ is central 

[□,div] = p,grad] = 0. 

In fact, □ is precisely the wave operator introduced quite some time ago by Lichnerowicz on the 
basis of its special algebra with gradient and divergence operators [9]. Finally, in the special case 
of constant curvature manifolds, choosing units in which the scalar curvature R = —n(n — 1), 
the curvature operator equals the sp(2, ]R) Casimir so that 

[div, grad] = □ - 2c. (8) 

If we include a further operator ord whose role is to count derivatives 

[ord,g] = [ord, N] = [ord,tr] = 0, 

[ord, grad] = grad, [ord, div] = div, [ord, □] = 2D, 

then, {ord, g, N, tr, grad, div, □} form a maximal parabolic subgroup of sp(4,lR) up to the 
rank 2 deformation by the sp(2, M) Casimir in (JSj) . On flat manifolds M, the operators {g, N, tr, 
grad, div, □} obey the Fourier-Jacobi Lie algebra of sp(2,]R). In the next Section, we present 
a novel reformulation of its universal enveloping algebra based on introducing a certain square 
root of the Casimir operator c. 

4 The Fourier— Jacobi algebra 

Let us first collect together the deformed Fourier-Jacobi Lie algebra built from geometric oper- 
ators on constant curvature spaces 





[tr,g 


] = 4N + 2n, 


[N,tr] 


= -2tr, 


[N,g] =2g, 


[N, grad] 


= grad 


[N, div] = -div, 


[tr, grad] 


= 2 div, 


[div,g] = 2 grad, 




[div, g 


rad] = □ - 2c. 



Its root diagram is given in Fig. [2j From now on we compute in the explicit realization given 
by its action on sections of the symmetric tensor bundle. Therefore we are working with linear 
operators, so any algebra we find is automatically consistent and associative. 
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N 



+2 g 



+1 



grad 



N,ord 



□ 



ord 



-1 



div 



-2 



tr 



+1 +2 

Figure 2. The root diagram for the Fourier- Jacobi Lie algebra. 



To start with we analyze the sp(2,R) Lie algebra built from {g, N,tr}. Unitary discrete 
series representations with respect to the adjoint involution tr' = g, = N are built from 
highest weights <3? such that 

N<J> = s<£>, tr$ = 0. 

The highest weight module is spanned by 

{$,g$,g 2 $,...}. 

We can characterize this representation by the eigenvalue s of N acting on the highest weight, 
or alternatively by the eigenvalue — s(s + n — 2) of the Casimir c acting on any state in the 
module. 

Conversely, given a eigenstate ^ of N and c, we can determine which discrete series repre- 
sentation it belongs to by repeatedly applying the trace operator 

tr fe ^ / = tr fc+1 ^, 

which implies that ^ can be expressed in terms of a highest weight vector as = g k <J>. Our key 
observation is that it is highly advantageous to introduce the linear operator k whose eigenvalue 
acting on $ is the depth k, namely 



M-C-l 



Hi 



(9) 



where 



n 



C 



n-2 



In other words, k measures how far the symmetric tensor ^ is from being trace-free. The 
operator 4AA is simply the right hand side of the tr, g commutator. More important is the 



8 



K. Hallowell and A. Waldron 



square root of the Casimir C which acts on the highest weight $ as 
n - 2" 



s + 



which explains equation 0. 

Our claim is that for many applications involving high powers of the operators g, N and tr, 
rather than normal or Weyl orderings it is far more expeditious to work with functions of Af, C 
(up to perhaps an overall power of g or tr). By way of translation, we note that a normal 
ordered product of g and tr can be expressed as PochhammeiH functions of (Af, C) : 

gtr = (Af-C- l)(N + C- 1), 

■( S tr) m -= S V = — fC ~ N+1 ^ fC+Af-1 

This claim has little significance until we introduce the doublet (grad,div). Indeed, since 
the Casimir c has a rather unpleasant commutation relation with either of these operators, 
computing in terms of (Af, C) may seem unwise. In fact this is not the case once one appropriately 
modifies the divergence and gradient operators. 

To motivate the claim we return to symmetric tensors. Suppose ip^ 1 ...^ La is trace-free, then 
its gradient V^^,,,,, s+1 ) is in general not trace free (unless the divergence of tp happens to 
vanish). Since we would like to work with states diagonalizing both c and N, it is propitious to 
replace the regular gradient with its trace-free counterpart 

We denote this operator by grad. Having introduced C and Af, it has the simple expression 

grad = grad — g div — r — , (10) 

Af + C-l v ; 

which we take to be its definition acting on any section of the symmetric tensor bundle. It is 
important to note that although this operator maps trace-free tensors to trace-free tensors, it 
is designed to maintain how far a more general tensor is from being trace-free. Therefore it 
does not project arbitrary tensors to trace-free ones. We also introduce a similar definition for 
a trace- free divergence following from the quantum mechanical adjoint grad^ = —div 

div = div — gradtr. (11) 

Af + C-l* v ; 

Note also that the linear operator Af + C — 1 is indeed invertible since its spectrum is 2s + 2k + 

n — 2 on eigenstates g fc< 3? (these expressions also make sense in dimensions n = 1,2 thanks to 

the operators div and gradtr). 

The beauty of the operators (div, grad) is that they commute with the depth operator k. 

This implies 

Afdiv = div (Af - 1), A/" grad = grad (Af + 1), 

CdTv = div(C - 1), C grad = grad (C + 1). 

Moreover, an easy computation using the definitions (|10[) and (jlip shows that the ordering 
of gradient and metric operators can be interchanged at the cost of only a rational function 
of (C,Af) 

~ Af + C-l ~ AA + C-1 ~^ 
grad g = g grad — -. , trdrv=- 7 divtr. 



3 



Recall that the Pochhammer symbol is defined as (x) m = x(x + 1) ■ ■ ■ (x + m — 1). 
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These relations allow us to invert equations (jlOp and ([IT 

, 1 M + C-3 • 1 M + C-3 ~ 

grad = grad H — - — — g div, 

s 2 C-l s 2C + C - 1 

l — AA + C-3 1 , AA + C-3 

div = - div — 1 — grad tr 



2 C-l 2 & C{N + C-l)' 

In turn we can now compute relations for reordering the gradient and trace operators 

tr grad = L_±£+iK_+^ 3 ) grad tr, 

~Tf ~ (AA + C + l)(AA + C-3) 

dlV g = g dlV — ; -7; . 

S S (N + C-l) 2 

The final relation we need is for div and grad. After some computations we find 

~~ C 2 {N + C + l)(N + C-3) 2 

div grad = — — — — — r grad div 

s (C + l)(C-l)(Af + C-l) 3 8 



+ 



2C 2 (AA + C + 1) 
(C + l)(Af + C-l) 2 



n \ n 



□ + 2[C + --lj (C-- + 1 



This result is valid in constant curvature spaces. The term in square brackets equals [□ — 2c] 
and will be modified accordingly upon departure from constant curvature. 

We denote the new algebra built from {g, M, tr, C, grad, div, □} by U(sp(2, M.) J ) and have 
collected together its defining relations in Fig. [3l As it is defined by linear operators acting 
on symmetric tensors, associativity is assured. An interesting, yet open, question is whether 
it can be defined on the universal enveloping algebra of sp(2,M) J . Nonetheless, the explicit 
symmetric tensor representation guarantees its consistency and therefore we may study it and 
its representations as an abstract algebra in its own right. 



5 Conclusions 

We have presented a detailed study of symmetric tensors on curved manifolds. The key technol- 
ogy employed is the quantum mechanics of a bosonic spinning particle model. Also, many of our 
constructions were originally motivated by studies of higher spin quantum field theories |10j . 
The spinning particle model presented here is one of a general class of orthosymplectic spin- 
ning particle models that describe spinors, differential forms, multiforms, and indeed the most 
general tensor-spinor fields on a Riemannian manifold [3]. 

There are many applications and further research avenues. One simple question is that 
given the strong analogy between the theory of differential forms and the symmetric tensor one 
presented here, are there symmetric tensor analogs of de Rham cohomology? The answer is yes. 
Recall, for example, the Maxwell detour complex 

d 5 
— c A°M — ► A 1 M -» A X M — ► A°M — ► 0. 



5d 

This is mathematical shorthand for Maxwell's electromagnetism in curved backgrounds. The 
physics translation is to replace the sequence of antisymmetric tensor bundles (A°M, A X M, A 1 ^^, 
A°M) by the words 

(gauge parameters, potentials, field equations, Bianchi identities). 
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CM = MC 
trg = (M + C + l)(M -C + l) 
g tr = (M + C-l)(M -C-l) 



Mtv = tr(M - 2) 
Ctr = trC 

Mdiv = div (M - 1) 
C div = div (C - 1) 

M + C-l 



gradg = g grad 



Mg = g(M + 2) 
Cg = gC 

M grad = grad (M + 1 
C grad = grad (C + l) 

M + C-l 



M + C + l 



tr div 



tr grad 



(M + C + l)(M + C-3) 



M + C + l 
grad tr 



div tr 



div g = g div 



(M + C-l) 2 
(M + C + l)(M + C-3) 
(M + C-l) 2 



C 2 (M + C + 1)(M + C-3) 2 — ^ ~ 
div grad = ln i 1W/ ^ tttt? ; /-> TTT g ra d div 



+ 



(C + 1)(C-1)(AA + C-1)3 
2C 2 (AA + C + 1) 
(C + 1)(AA + C-1) 



□ + 2[C + - 



C -2 +1 



Figure 3. Defining relations for the W(sp(2,]R) J ) algebra. 



Then the fact that (|12p is a complex implies that Maxwell's equations 5<i^4 = are gauge 
invariant because Sdd = 0, and subject to a Bianchi identity as 55d = 0. An analogous complex 
exists for symmetric tensors although not in general backgrounds, for brevity we give the flat 
space result [12l [10] 



grad 

SM — ► SM 



div 

SM — ► SM 

T 



o. 



G 



where 



G 



□ — grad div H — (grad 2 tr + g div 2 ) g ( □ H — grad div ) tr. 



Notice that we have specified no grading on the symmetric tensor bundle SM. In fact the 
operator G is the generating function for the equations of motion (and actions) for massless 
higher spins of arbitrary degree. A very fascinating question is whether such complexes exist 
for the most general orthosymplectic spinning particle models - preliminary studies suggest an 
affirmative answer |16j . 

Another interesting open question is the generality of the algebra IA (sp(2, R) J ). For example, 
does there exist an algebra U (osp(2p\Q) J ) where osp denotes the orthosymplectic superalgebra. 
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Since the key idea is to include the square root of the Casimir operator in the algebra, higher 
rank generalizations ought involve the higher order Casimir operators. A positive answer to this 
question would be most welcome and is under investigation [17| . 
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